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We address the issue of the dynamical origin of scale-free link distributions. We study 
a two-dimensional lattice of cooperatively interacting units. Although the units interact 
only with the four nearest neighbors, a sufficiently large cooperation strength generates 
dynamically a scale-free network with the power law index v approaching 1. We explain 
this result by using a new definition of network efficiency determined by the Euclidean 
distance between correlated units. According to this definition the link distribution 
favoring long-range connections makes efficiency increase. We embed an ad hoc scale-free 
network with power index v > 1 into a Euclidean two-dimensional space and show that 
the network efficiency becomes maximal as v approaches 1. We therefore conclude that 
v = 1 emerging from the cooperative interaction of units may be a consequence of the 
principle of network maximal efficiency. 

© 2011 Elsevier B.V. All rights reserved. 


1. Introduction 

The work of Albert and Barabasi (AB) [1] prompted interest in scale-free networks, namely networks with the inverse 
power law distribution of links 

P(k) oc 1, (1) 

with k denoting the number of links per node. The AB model assumed the preferential attachment of the Yule process and 
gave rise to the power-law index v = 3. The ubiquity of these networks has led to the research into the dynamical origin 
of the scale-free property. Fraiman et al. [2] has recently studied the two-dimensional Ising model at criticality and using 
a suitable threshold p + established a link between the nodes with a dynamical correlation exceeding that threshold. The 
interesting result of that procedure was the emergence of the scale-free distribution density of Eq. (1) with a power index 
v that, with the proper choice of average connectivity ( k >, gives values close to v = 2 [3], It is remarkable that the same 
procedure, applied to the brain, yields results very close to those afforded by the Ising model [2], 

Like the work of Ref. [2], we begin with a regular two-dimensional network, but then take an entirely different approach. 
Rather than using the Ising model, we adopt the Decision Making (DM) model of Refs. [4,5] for the creation of the dynamically 
induced network. The key result is that while we also find the emergence of a scale-free network at criticality, our network 
has an exponent v close to 1. 
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This finding suggests that the cooperative interaction between the units in the network generated from the DM model 
should make the network very efficient. However, we find herein that the conventional measures of network efficiency [6] 
are relatively insensitive to v < 2. The second result of this paper is a new measure of network efficiency, the perception 
length, that overcomes this limitation, and shows that network efficiency increases as v approaches 1. 

The outline of this paper is as follows. We devote Section 2 to proving the dynamical emergence of the scale-free network 
with v = 1. Section 3 shows that v approaching 1 maximizes the perception length while realizing a statistical condition 
equivalent to the dynamically generated network. We devote Section 4 to concluding remarks. 


2. Dynamically induced scale-free network with v = 1 


The question of whether the network emerging from dynamical correlations is scale-free is not trivial. For this paper, we 
adopt the Decision Making (DM) model of Refs. [4,5]. Although the DM model yields a phase transition very similar to that 
of the conventional Ising model [7], the DM model does not have a Hamiltonian foundation. Furthermore, while the Ising 
model rests on the action of a thermal bath at a finite temperature, and the single spins in isolation and with no thermal 
driving would be dynamically frozen, the units of the DM model, in isolation, are driven by Poisson dynamics. This difference 
may lead to a dynamically created scale-free network different from that of Ref. [2], 

The DM model considers N discrete variables located at the nodes of a two-dimensional square lattice. Each unit s, is a 
stochastic oscillator and can be found in either of two states, “yes” (+1) or “no" (—1). The cooperation among the units is 
realized by setting the transition rates between two states to the time-dependent form: 


gu(t) =ge K 


M 2 (f)-M 1 (t) 
M 


( 2 ) 


and 


g 2 i(t) =ge K 


M 2 (t)-M 1 (t) 

M 


(3) 


where M denotes the total number of nearest neighbors, which in the case of a two-dimensional lattice results in M = 4. 
Mi(t) and M 2 (t) are the nearest neighbors who are making the decision “yes” and “no”, respectively. The single unit in 
isolation, I< = 0, fluctuates between “yes" and “no”, with the rate g. When coupling constant I< > 0, a unit in the state “yes” 
(“no") makes a transition to the state “no” ("yes”) faster or slower according to whether M 2 > Mi (Mi > M 2 ) or M 2 < Mi 
(Mi < M 2 ), respectively. The quantity K c is the critical value of the control parameter I(, at which point a phase-transition 
to a global majority state occurs. It can be shown that for a lattice of infinite size I( c = 2 ln(l + V2) ~ 1.7627. 

All simulations are implemented on a lattice of JV = 100 x 100 nodes with periodic boundary conditions. In a single time 
step, a run over the entire lattice is performed and for every node an appropriate transition rate (g J2 org 21 ) is calculated 
according to which a node is given possibility to change its state. Under those conditions, setting parameter g = 0.01, the 
phase transition to the global majority case occurs at K c ~ 1.70. To study the dynamically induced network topology, we 
considered the DM model with this critical value of coupling constant. After an initial 10 6 time steps, we record 2000 lattice 
configurations, obtaining the dynamics of each node {s, (t)}. In the next step, we evaluate the linear correlation coefficient 
between the i-th and the j-th node: 


,. .. (Si(t)Sj(t)) - (Si(t)) (Sj(t)) 

r(ij) = , - (4) 

yj (S?(0) - (S,(t)) 2 v /(S J 2 (t)) - <Sj(t)) 2 

where (• • •) stands for the time average. Nodes i and j of the dynamically induced network are connected by a link when the 
correlation between nodes i and j of the two-dimensional lattice is greater than a given positive threshold p + . 

We investigated a wide range of positive thresholds p + = [0.10 : 0.90], As expected, correlation networks constructed 
with increasing value of the threshold p + are characterized by a decreasing number of links, which results in significant 
changes in the corresponding degree distribution. Networks obtained with low values of the threshold p + result in narrow, 
peaked degree distributions and relatively high values of the mean degree {!<). Topologies resulting from very strong 
correlations, however, are characterized by a distribution of links that has a mean close to the minimum value of (k) = 1 
and drops off very rapidly. Interestingly, in between those conditions there exists a correlation threshold that creates a 
dynamically induced network with a power-law distribution of degrees that is shown in Fig. 1. The connectivity of a network 
obtained under this particular condition is inspected and only nodes that are part of a giant component are considered to be 
a part of new topology. The resulting network consists of 6740 nodes with a mean degree of (k) = 38. For a correlation 
threshold p + = 0.61, we find that the distribution of links follows an inverse power law with the scaling parameter 
v = 1.20 ± 0.23 (R 2 = 0.9556). 

How can the dynamical origin of this extremal condition be explained? From an intuitive point of view there exists a 
connection between the scale-free distribution of links and Zipfs law [8,9]. In fact, if we rank the nodes moving from the 
richest to the poorest, according to Zipf [10] we should obtain the relation between the connectivity k and rank r: 


C 

r a ’ 


k(r) = 


(5) 
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Fig. 1. Degree distribution p(k) for correlation network obtained with a positive threshold p+ = 0.61. The line follows estimated power-law relation with 
the exponent of v = 1.20. 

with the power-law index a being very close to 1. If we select randomly a given rank r with probability ;r(r)dr of being in 
the interval [r, r + dr], this should be equal to the probability p(/<)d/<, of having fc links in the interval [k, k + dfc], 

p(k)dk = ;r(r)dr. (6) 

Assuming that tz(t) is constant, we obtain after some algebra 

p(k) oc — 1 - r , (7) 

k 1+ ? 

and thus the two power-law indices are related by 

v = l + -. (8) 

a 

This leads us to conclude how an efficient complex network is organized if Zipfs law applies. For example, this suggests 
that the intelligence of human brain should be characterized by v = 2. The recent research work of Fraiman et al. [2] seems 
to confirm this prediction. The authors of Ref. [2] studied the two-dimensional Ising model at criticality and using a suitable 
threshold p + established a link between the nodes with a dynamical correlation exceeding that threshold. The interesting 
result of that procedure was the emergence of the scale-free distribution density of Eq. (1) with a power index v that, with 
the proper choice of average connectivity (fc), gives values close to v = 2 [3], It is remarkable that the same procedure, 
applied to the brain, yields results very close to those afforded by the Ising model [2], 

It is also interesting to notice that Boettcher and Percus [11] applied the method of extremal optimization to a ising- 
like model to evaluate the rank distribution of Eq. (5) and found that when the number of units tend to infinity the rank 
distribution tends to a = 1. This result implies on the basis of the earlier remarks that the cooperation of units in an Ising- 
like model is expected to confirm Zipfs prediction and consequently v = 2. For results confirming this prediction, see also 
Refs. [12,13], 

We have found, therefore, that the DM model generates a scale-free network at criticality, but with v approaching 1 
which violates Zipfs condition of a = 1. The emergence of a 1 corresponds to the number of links k dropping very 
quickly when we move from r = 1 to subsequent ranks. 

3. Embedding scale-free networks in a two-dimensional regular lattice 

The efficiency of a network has been mainly studied on the basis of topological length L [6]. The topological distance 
between two nodes is the minimum number of steps needed to move from one node to another. The quantity L is the 
corresponding mean value of the minimum topological distance between nodes. As a consequence, the emergence of scale- 
free networks can be explained by making the conjecture that real networks evolve in time so as to realize maximal 
efficiency. In the case of the scale-free network this efficiency is measured by the minimal value of L. 

In 2003, Cohen and Flavlin [14] determined that scale-free networks are very efficient. They calculated that networks 
with N nodes have a topological length L given by 

L ~ ln(lnN), (9) 

for a power-law index v < 3, which is smaller than in the AB model. 

As a consequence of Eq. (9), the topological length L remains small even if N is very large. This relation between L and N 
suggests that the complex networks with v < 3 are very efficient. 

The region v < 2 has not been studied by many authors. As an interesting example of earlier work in this region we 
refer to the work of Refs. [15,16], In both cases the authors adopt a kind of generalization of the perspective of Albert and 
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Fig. 2. Perception length P as a function of v for ad hoc scale-free networks. The circles represent networks with constant (k) = 4.95 ; the squares represent 
networks with constant N = 1313. 

Barabasi [1], namely, they grow a network with special prescriptions allowing them to overcome the limiting condition 
v > 3 and yielding v < 2. 

More recently there has been some interest in the efficiency of hierarchical networks with v < 2[17]and v ~ 1 [18,19], 
These networks are characterized by very small topological lengths. In the case of Ref. [19], forN —> oc the authors predict 
the value L = 2. Our numerical calculations show that L is insensitive to v in the range 1 < v < 2. Thus, there are good 
reasons to believe that the prediction of Eq. (9) applies with no significant dependence on v. 

There are also good reasons to believe that the dynamic source of scale-free networks has a probabilistic rather than 
a deterministic origin, as we shall demonstrate. The work of Fraiman et al. [2], although fitting Zipfs condition, allows us 
to define the concept of perception length, which is then used to shed light onto the condition v = 1. The nodes of a 
two-dimensional regular lattice, where each unit interacts with only its four nearest neighbors, correspond to a network 
whose topological efficiency is very low. However, if we let the units cooperatively interact, and we set a link between two 
nodes, then when their dynamical correlation is large enough, the resulting dynamical network may have a much lower 
L. We record the Euclidean distance between each pair of correlated nodes, and we define the network efficiency as the 
corresponding mean value. More precisely, we define the network efficiency by means of the quantity 



( 10 ) 


where 



( 11 ) 


and dy is the Euclidean distance between node i and its nearest neighbor j. We refer to P as the perception length of the 
system, this being the measure of the network efficiency replacing 1 /L For numerical analysis, the perception length P is 
calculated as follows. 

1. Create a 2-dimensional square lattice, where a is the length (in units of the lattice) of the side of the square, and it is 
necessary that >/n < a. 

2. Embed the network on the lattice by randomly assigning the network nodes to lattice points. 

3. Using periodic boundary conditions for the lattice, calculate the Euclidean length of each link. 

4. These lengths are then used to calculate P in accordance with Eq. (10). 

This procedure was then applied to scale-free networks generated according to the method of Catanzaro et al. [20] with 
(k) = 4.95 and values of v ranging from 1 to 2. The nodes were embedded randomly in a square lattice of size just sufficient 
to contain the network with the largest number of nodes. The perception lengths were calculated, and the results show that 
over this range the perception length increases as v approaches 1 (Fig. 2). 

We now compare the sensitivity of the conventional network efficiency (defined as the reciprocal of the network length 
1/L) and the perception length P to the variation in v. These results are shown in Fig. 3. Note that the perception length P and 
1 /L have been normalized to allow the presentation in a convenient scale. These results show that the perception length, as 
a measure of network efficiency, is more sensitive to changes in v than the conventional 1 /L 

Applying this approach to the DM network discussed in Section 2 above yields a mean perception length = 1382. To 
provide a basis for comparison, we calculated the perception length for a comparably sized ad hoc scale-free network, also 
prepared in accordance with the algorithm in Ref. [20]. The results are compared in Table 1, which shows a very good 
agreement. 
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Fig. 3. Network efficiency as a function of v. In both figures, the circles represent the normalized perception length, and the squares represent the 
normalized 1/L Error bars are not shown; standard deviations for the normalized perception lengths were very small (of the order of 0.005). 


Table 1 

Comparison of perception lengths for the DM network and a comparable 
ad hoc scale-free network. 



DM network 

Ad Hoc network 

Degree distribution exponent, v 

1.20 

1.19 

Number of nodes, N 

6740 

6744 

Number of links, k 

128,448 

128,413 

Perception length, P 

1382 

1373 


4. Concluding remarks 

The main result of this paper is the proof that the perception length P defined by Eq. (10) is more sensitive to v in the range 
1 < v < 2, and may therefore be a more useful measure of network efficiency than the topological length L [6], This new 
definition is based on the assumption that the scale-free networks, even those derived using deterministic and hierarchical 
arguments [17-19], are compatible with the dynamical derivation. The apparent conflict between the two kinds of network 
is explained by noticing that the correlation between different nodes evaluated according to the numerical prescription 
of Section 2 rests on time windows of finite size and moving these windows to different time regions generate different 
networks with the same hierarchical topology. An alternative but probably more attractive way to explain this is that the 
leadership moves in time from some nodes to others. The nodes of the two-dimensional regular lattice that we used for the 
dynamical derivation of the free-scale distribution of Eq. (1) are totally equivalent, and the rank differences between the 
richest and the poorest nodes are temporary and in the long time run the same high rank is shared by all the nodes. 

There are some questionable aspects of the dynamical derivation of v = 1, realized in Section 2. First of all this result 
depends on an arbitrarily selected threshold p + . In addition, the power index v derived from the data in Fig. 1 is affected 
by the large fluctuations of the distribution density p(k) for large values of k, thereby explaining why we derive v = 1 
rather than v = 2, as in the work of Ref. [3], However, on one hand, this observation led us to the discovery of the new 
measure of network efficiency of Eq. (10) and, on the other hand, as already mentioned in Section 1, the DM model is not 
quite equivalent to the Ising model. 

Furthermore, the interesting fact that a = 1 (Zipfs law) is related to v = 2 is not totally convincing. The authors of 
Ref. [9] established a connection between topological and dynamical properties by running a random walker on a scale-free 
network with the power index v. They found that the waiting time distribution density for returns of the walker to a given 
node are described by an inverse power law with power index /z = 3— v. Using a generalized central limit theorem, they also 
proved that the same index /z also describes the inverse power law form of the probability density of relative frequencies 
p(f) ~ / -M , so that, in analogy with Eq. (8), a rank-frequency law/ ~ r _ “ yields /z = 1 + 1/a, hence a = 1/(2 — v), 
namely a = 1 for v = 1. It is interesting to note that v = 1 yields /z = 2, and that the value /z sa 2 emerges from the recent 
neurophysiological work of Refs. [21,22]. 

Thus, in addition to the new efficiency measure of Eq. (10), which is the unquestionable result of this paper, we are also 
led to make the plausible conjecture that the DM model may be an adequate tool to study the emergence of intelligence 
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in the sense pointed out by Couzin [23] and Cavagna et al. [24]. According to these authors the flocks of birds, as well as 
the brain [25,26], are systems at criticality, with a correlation length as large as the size of these systems. We are also led 
to make the conjecture that the dynamical generation of a scale-free network with v = 1 may be the manifestation of the 
principle of network maximal efficiency: a network of cooperatively interacting units evolves naturally so as to maximize 
P, the new measure of network efficiency. 

References 

[1] R. Albert, A.-L. Barabasi, Statistical mechanics of complex networks, Rev. Modern Phys. 74 (2002) 47-97. 

[2] D. Fraiman, P. Balenzuela, J. Foss, D.R. Chialvo, Ising-like dynamics in large-scale functional brain networks, Phys. Rev. E 79 (1-10) (2009) 061922. 

[3] E. Tagliazucchi, D.R. Chialvo, The collective brain, in: P. Grigolini, B.J. West (Eds.), Decision Making: A Psychophysics Application of Network Science, 
World Scientific, Singapore, 2011. 

[4] S. Bianco, E. Geneston, P. Grigolini, M. Ignaccolo, Renewal aging as emerging property of phase synchronization, Physica A 387 (2008) 1387-1392. 

[5] M. Turalska, M. Lukovic, B.J. West, P. Grigolini, Complexity and synchronization, Phys. Rev. E 80 (1-12) (2009) 021110. 

[6] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, D.U. Hwang, Complex networks: structure and dynamics, Phys. Rep. 424 (2006) 175-308. 

[7] M. Turalska, B.J. West, P. Grigolini, Temporal complexity of the order parameter at the phase transition, Phys. Rev. E 83 (1-6) (2011) 061142. 

[8] L.A. Adamic, B.A. Huberman, Zipfs law and the Internet, Glottometrics 3 (2002) 143-150. 

[9] P. Allegrini, P. Grigolini, L. Palatella, Intermittency and scale-free networks: a dynamical model for human language complexity, Chaos Solitons Fractals 
20(2004)95-105. 

[10] G.K. Zipf, Selected Studies of the Principle of Relative Frequency in Language, Harvard University Press, Cambridge, MA, 1932; Human Behavior and 
the Principle of Least Effort, Addison-Wesley, Cambridge, MA, 1949. 

[11] S. Boettcher, A.G. Percus, Optimization with extremal dynamics, Phys. Rev. Lett. 86 (2001) 5211-5214. 

j 12] S. Boettcher, Acquiring long-range memory through adaptive avalanches, in: P. Grigolini, B.J. West (Eds.), Decision Making: A Psychophysics 
Application of Network Science, World Scientific, Singapore, 2011. 

[13] It is important to notice that in the work of Ref. [12] the role of Eq. (5) is played by Eq. (4) that refers to the probability that a node of rank r is updated. 
However, if we divide Eq. (5) by the maximum number of links that a node may have, k is replaced by quantity that can be interpreted as a probability 

[14] R. Cohen, S. Havlin, Scale-free networks are ultrasmall, Phys. Rev. Lett. 90 (1-4) (2003) 058701. 

j 15 ] H. Zhou, Scaling exponent and clustering coefficients of a growing random network, Phys. Rev. E 66 (1-6) (2002) 016125. 
j 16] H. Seyed-allaei, G. Bianconi, M. Marsili, Scale-free networks with an exponent less than two, Phys. Rev. 73 (1-5) (2006) 046113. 

[17] A.-L. Barabasi, E. Ravasz, T. Vicsek, Deterministic scale-free networks, Physica A 299 (2001) 559-564. 

j 18] M. Chen, B. Yu, P. Xu, J. Chen, A new deterministic complex network model with hierarchical structure, Physica A 385 (2007) 707-717. 

[19] Y. Gua, J. Suna, A tree-like complex network model, Physica A 389 (2010) 171-178. 

[20] M. Catanzaro, M. Boguna, R. Pastor-Satorras, Generation of uncorrelated random scale-free networks, Phys. Rev. E 71 (1-4) (2005) 027103. 

[21] P. Allegrini, D. Menicucci, R. Bedini, L. Fronzoni, A. Gemignani, P. Grigolini, B.J. West, P. Paradisi, Spontaneous brain activity as a source of ideal \/f 
noise, Phys. Rev. E 80 (1-13) (2009) 061914. 

[22] P. Allegrini, D. Menicucci, R. Bedini, A. Gemignani, P. Paradisi, Complex intermittency blurred by noise: theory and application to neural dynamics, 
Phys. Rev. E 82 (1-4) (2010) 015103 (R). 

[23] I. Couzin, Collective minds, Nature 445 (2007) 715. 

[24] A. Cavagna, A. Cimarelli, I. Giardina, G. Parisi, R. Santagati, F. Stefanini, M. Viale, Scale-free correlations in starling flocks, Proc. Natl. Acad. Sci. 107 
(2010)11865-11870. 

[25] D.R. Chialvo, Emergent complexity: what uphill analysis or downhill invention cannot do, New Ideas in Psychology 26 (2008) 158-173. 

[26] D.R. Chialvo, Complex emergent neural dynamics, Nature Physics 6 (2010) 744-750. 


